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D script! n 

La presents invention a pour objet un proc6d6 d f excitation par un signal de radiofrequence dans une expe- 
rimentation de resonance magnetique nucl£aire (RMN). Elie trouve plus particullferement son application dans 

s !e domaine medical ou des machines de RMN sont utilises a des fins de diagnostic chez I'homme. Elle peut 
neanmoins trouver son application dans les domaines du laboratoire ou de I'industrie des que des ph£nom&nes 
de resonance magnetique sont mis en jeu. 

Le phenomene de la resonance magnetique nucleate resulte de I'orientation prise par les moments 
magnetiques des spins nucleates, ci-apres denomm6s particules, d'un corps quand ce corps est soumis a un 

10 champ magnetique orientateur constant. Pour faire apparaitre ce phenomene on excite le corps par un signal 
de radiofrequence. On peut mesurer le signal de desexcitation des particules quand elles retoument a ieur £quf- 
libre, quand Ieur moment magnetique se reoriente avec le champ orientateur, a la fin du basculement de cette 
orientation dependant de I'excitation. L'excitation est appl iqu6e a la totality du corps soumis a examen , et toutes 
les particules du corps 6mettent un signal de desexcitation a son issue. Dans le domaine medical ou on cherche 

15 a representor des images de coupes du corps examine (d'un patient), on a coutume de cantonner la naissance 
du ph6nomene de resonance & des tranches s6lectionn6es dont on veut donner Pimage. Dans ce but on modifie 
localement les conditions de resonance de manure a ce que seule ia tranche s£iectionn6e se trouve dans des 
conditions de resonance correspondant & la frequence de l'excitation radiofrequence. En dehors de la tranche, 
les moments magnetiques ne sont pas bascules, les particules exterieures a la tranche ne restituent aucun 

20 signal de desexcitation. 

La modification des conditions locales de resonance s'obtient en augmentant, ou en diminuant, I'intensite 
du champ orientateur dans la tranche s6iectionn6e. Pour des raisons scient'rfiques et technologiques il n'est 
pas possible de modifier brutalement I'intensite d'un champ magnetique continu de part et d'autre d'une tranche 
de I'espace. La modification du champ prend alors la forme d'une variation reguliere en fonction de I'abscisse 

25 des lieux de I'espace sur un axe perpendiculaire a la tranche £ seiectionner. L'intensite du champ orientateur 
subit doncseton cet axe un gradient, etpratiquement les champs suppiementaires, additionnes au champ cons- 
tant et homogene dans tout I'espace pour provoquer cette variation, sont commun6ment ap pell 6s les gradients 
de champ. La frequence de resonance du phenomene de resonance est proportionnelle a I'intensite du champ 
orientateur. Pour seiectionner une tranche il convient done, dans l'excitation radiofrequence, que le spectre 

30 de l'excitation ne contienne des composantes que comprises dans une bande lirnitee par des frequences de 
coupures correspondant aux frequences de resonance des particules des bonds de la tranche. 

L'obtention d'un signal d'excitation a spectre etroit entraine que ce signal ne peut 6tre instantane. II dure 
un certain temps. En pratique des puissances de quelques kilowatts, appliquees pendant des durees de quel- 
ques dizaines de millisecondes, sont courantes. Le fait que I'impulsion d'excitation n'est pas instantan6e et 

35 qu'au contraire elle dure un certain temps provoque alors des perturbations dans le signal de desexcitation. 
En etTet, du fait du gradient de champ, les particules des bords de la tranche ne resonnent pas a des m§mes 
frequences. Par exemple si le milieu de la tranche resonne a une frequence f 0r les bords amont et aval de la 
tranche, par rapport £ I'orientation du champ, peuvent resonner a des frequences f 0 - Af. ou f 0 + Af respects 
vement et il en r6sulte que I'historique de l'excitation dans les bords de la tranche n'est pas le meme qu'au 

40 centre bien que la dur6e de l'excitation ait 6te la meme pour toutes les particules de la tranche. Si le bascule- 
ment geometrique de tous les moments magnetiques est conrectement obtenu pour une 6nergie d'excitation 
donn6e (par exemple ils se retrouvent tous dans un plan perpendiculaire au champ orientateur pour une exci- 
tation a 90°), Ieur orientation n'est alors pas unique. Elle peut au contraire fitre distribuee dans toutes les direc- 
tions contenues dans un plan. L'orientation de chaque moment y depend de I'abscisse des points considers 

45 dans la tranche. En effet dans Ieur mouvement de precession forcee pour rejoindre le basculement geometrique 
recherche, les moments magnetiques ont emprunte des trajectoires differentes. Ces trajectoires sont diff6ren- 
tes car les vitesses de precession sont Ii6es aux frequences de resonance. Comme I'orientation initiale des 
moments magnetiques est la meme pour tous, on comprend qu'il n'y a aucune raison a priori pour que cette 
orientation soit unique a la fin du basculement. Si on n'y prend garde les moments magnetiques des divers 

50 points de la tranche peuvent se trouver a ce point desorientes, dephases, les uns par rapport aux autres que 
le signal eiectromagnetique de restitution peut avoir une moyenne nulle. II est alors impossible d'en tirer une 
information physique. 

En butre I ph6n mfcne d resonance magnetique ne s situ pas a une frequ nee uniqu , mais plutfit a 
des frequenc s voisines. L'ecart entre ces frequences et un frequ nee de reference est appeie le decalage 
55 chimlque (en literature anglo-saxonne : chemical shift). Cet ecart est en general inferieur a 2 Af. II est recherche 
par les praticiens une excitation qui seraitcapabl ,p urune xp6rimentationdonn£e,dedifferencierlar£p nse 
correspondant a chaque valeur de se decalage chimlque. De cette maniere il serait possible, en particuli r dans 
les images medicales, de faire une separation eau-graisse dans les tissus. En ffet pour les diff£rents tissus, 
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si le d6calage de I'eau par rapport a la graisse a una meme valeur n frequence, lee amplitudes des compo- 
santes n* nt pas des memes va! urs d'un tissus a I'autr . Par ailleurs, dans un autre type d'applicati n t n 
particulier dans le domalne d I'imagerie, il est envisage de construire des excitations dites multicoupes ou 
plutdt multitranches. L but s rait de profiler d'un gradient qui s'etend dans I'espace sur une I ngu ur bien 

5 supen'eure a Tepaisseur de la tranche et de disposer d'une excitation susceptible cf exciter, en meme temps, 
plusieurs tranches, non contigues de preference, de maniere a accelerer les procedes de reconstruction 
d'image. Or le travail entrepris pour la definition des impulsions d'excitation classique (a 90° et a 1 80°) a paru 
a ce point complexe que I'homme du metier n'envisage pas d'en augmenter encore la complexite pour produire 
maintenant de telles excitations. r 

10 Dans I'etat de la technique le probleme de la determination du signal d'excitation a ete resoiu par des 
approximations successives et par des tatonnements qui peu a peu ont permis de determiner Thistorique de 
son amplitude pendant la dur6e de impulsion de maniere a ce que des resultats desir6s soient atteints. Cepen- 
dant, malgre les mois de travail necessaires pour arriver a ces resultats, {'optimisation des impulsions d'exci- 
tation n'est pas garantie, notamment en ce qui concerne la puissance totale mise en oeuvre ainsi que la 

15 puissance instantan6e. Le but de la presente invention est de proposer un procede de determination de I'exci- 
tation qui permette en plus de tenir compte d'un critere d'optimlsation , meme 6ventuellement different de ceux- 
ci. 

La technique utilisee dans I 'etat de la technique pour determiner les excitations en RMN de protons consists 
essentiellement a diviser une dur£e pr6determin6e d'excitation en un nombre suffisant de periodes, par exenv 

20 pie 40, et a calculer pour chacune des abscisses de la tranche la rotation, pendant chaque periode, du moment 
magnetique des protons qui s'y trouvent. Le caicul est entrepris en fonction d'un parametre representatif de 
Pamplitude de I'excitation pendant la periode correspondante. On fait evoluer dans le temps ce parametre pen- 
dant la duree de I'excitation. Le mouvement du moment magnetique est un mouvement de precession du a la 
frequence de resonance, done dependant de l'abscisse dans la tranche, combine a un mouvement de nutation 

25 impose par I'excitation radiofrequence, identique pour toute la tranche. Aussi il met en jeu pour determiner la 
composante selon une direction de I'espace du moment magnetique d'une particule a la fin d'une periode, les 
composantes selon les trois directions de I'espace du moment magnetique de cette particule avant le debut 
de la periode, ces trois composants dependant, du fait du gradient, de l'abscisse du point dans la tranche. En 
consequence le caicul doit dtre entrepris abscisse par abscisse pour toutes les abscisses dans I'epaisseur le 

30 long de la tranche. 

Dans la pratique ces calculs ont 6te entrepris pour une centalne d'abscisses reparties regulierement au 
travers, et un peu de part et d'autre, de ia tranche. Le travail de determination a consists a fairs la totalite de 
ce caicul, et a le re prendre un certain nombre de fois en modifiant, p6riode par periode, I'historique de revolution 
du parametre de maniere a obtenir en fin d'impulsion une definition recherchee en amplitude et en phase du 

35 basculement des moments magn6tiques des particules soumises a ('experimentation. Ces calculs concement 
toutes les particules dans la tranche, et voisine de la tranche, e'est a dire un ensemble de frequences de reso- 
nance. Dans une autre approche, justifiee uniquement pour de faibles angles de nutation, par exemple infe- 
rieurs a 30°, on utilise le fait que 1'expresslon de I'aimantation transverse en fonction de la frequence de 
resonance est directement la transformee de Fourier du signal tempore! d'excitation. Cette hypothese devient 

40 peu vaiide quand Tangle de nutation recherche vaut 90°. Elle est completementfausse quand il vaut 180°. 

La presente invention permet de rem6dier aux inconvenients cites en proposant une methods de caicul 
de ('impulsion d'excitation completement differente. Dans cette methode on ne calcule pas en chaque abscisse 
directement Tangle de nutation de chaque moment magnetique. On exprime seulement, sous une forme ana- 
lytique polyn6miale en fonction de la frequence de resonance, la probability que ce moment ait une orientation 

45 desiree a la fin de I'excitation. L'expression sous forme polynomial permet de s'affranchir, dans les calculs, 
de Tinformation relative a la localisation des protons concemes (ou de leur frequence de resonance ce qui 
revient au mSme) de maniere a pouvoir proposer une solution presentant pour toutes les frequences de reso- 
nance impiiquees un meilieur compromis vis a vis d'un critere d'optimisation. En ramenant cette optimisation 
a des syntheses de fDtres par transformee de Fourier, Tinvention permet alors d'aboutir rapidement a une solu- 

so tion recherchee (en quelques minutes). 

La presente invention a pour objet un procede d'excitation radiofrequence dans une experimentation de 
RMN, excitation en reponse a laquelle les bascuiements de moments magn6tiques de particules soumises a 
('experimentation sont d6finis en amplitude et en phase pour un ensemble de frequences de resonance donn6 
caracterise en ce que : 

55 — necritl'equationdeSchrddingerr lative a I'amplitud d probalitedelafoncti nd' nd d spins reson- 
nant a un frequence de. resonance donn6e et dont Thamiltonien depend de la frequence et de I'excitation 
radiofrequence qu'ils subissent, 
— on int6gre cette equation par rapport au temps, 
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— on en deduit le xpressions analytiques polyndmial s des amplitudes d probability k la fin de P xci- 
tation en fonction de !a frequence de resonance, ces expressions po!yn6miales etant de degr6 6gal k un 
nombre de period s 6l6mentaires dont la somme des durees est egale k la duree de Pexcitation, et ces 
expressions etant parametrees par des coefficients repr^sentatifs de Phistorique de Pexcitation pendant 

5 ces periodes, 

— on optimise ces coefficients k Paide d'algorithmes de calcul de facon k ce que les polyndmes s'appro- 
chent au mieux d'une forme desir6e 

— on d6duit les caracttrisitiques, pendant chaque pSriode, de Pexcitation qui conduit k cette r£ponse defi- 

nle. 

10 — et on excite les particules avec une telle excitation. 

La pr£sente invention sera mieux comprise k la lecture de la description qui suit et k Pexamen des figures 
qui Paccompagnent. Cette description n'est pas limitative de ('invention. Notamment, Papplication d'imagerie 
6voqu6e n'est qu'un cas particulier de determination d'une excitation radiofrequence. 
Les figures montrent : 

15 — figure 1 : une representation schematique d'une machine de RMN ; 

— figures 2 et 3 : des allures d'enveloppes d'amplitude d'excitation d6termin6es selon Pinvenb'on. 
La figure 1 repr6sente schematiquement une machine d'imagerie par resonance magnetique nucleaire. 
Cette machine comporte essentiellement des moyens 1 pour produire un champ magnetique orientateur, cons- 
tant et homogene, B 0 dans un lieu ou est place un corps 2 k examiner. Le corps 2 est soumis k cet endroit k 

20 une excitation par un signal de radiofrequence, produite par !e generateur 3 alimentant une antenne 4. 
L'antenne sert egalement k preiever le signal de dgsexcitation et k Pacheminer, au moyen d'un dupfexeur 5, 
sur des moyens de reception et de traitement 6. A Pissue du traitement, Pimage d'une coupe 7 du corps 2 peut 
etre visuallsee sur un dispositif 8 de visualisation. Pour soumettre les moments magnetiques des protons du 
corps 2 situ£s dans la tranche 7 au phenomena de resonance, et eux seuls, on a mis en oeuvre des bobines 

25 de gradients 9 alimentees par un generateur 10 d'impulsions de gradient Tous ces moyens fonctionnent sous 
le pilotage d'un s6quenceur 11. En supposant que le champ magnetique sup!6mentaire applique par les bobi- 
nes 9 varie li n6airement selon un axe z parallfele au champ B 0 . on peut dire que les bords de la tranche 7, bords 
12 et 13 respectivement en amont et en aval, rencontrent des conditions de resonance differentes et que les 
frequences de resonance des protons qui s'y trouvent sont differentes. 

30 Dans i'inveniion, pour determiner Pexcitation radiofrequence, on ecrit Pgquaticn ds Schrcdinger relative k 
Pamplitude de probabtlite de la fonction d'onde de spins resonnant k une frequence de resonance donnee. 
Cette expression est la suivante : 

Dans cette formule H a est Phamiltonien de spin, il depend de la frequence et de Pexcitation radiofrequence 
que subissent les spins r6sonnants. Le terme | <p > repr6sente ia fonction d'onde. Cet hamiltonien peut s'6crire 
40 de la maniere suivante : 



45 



H A 
a 2 



x ' y 



x ' y 



Dans ces equations j est ie nombre comptexe tel que j 2 = - 1. Les pulsations w, et w x et w y sont, respec- 
tivement, I'ecart de pulsation de resonance dans ia tranche 7 en fonction de Pabscisse selon z des protons 
so concernes, et les representations des composantes radiofrequence du champ d'excitation Bj 6mis par 
Pantenne 4. Ces composantes sont evaiuees dans un repfere tournant k une pulsation w 0 , correspondant k la 
porteuse de Pexcitation et k la frequence de resonance des protons situes au centre de la tranche 7. De Pequa- 
tion de SchrOdinger on tire I' xpression survant : 

55 
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10 Oil a et b sont !es ampl itudps de probabilite pour que les spins des protons considered soient dans un etat 
respectivement positif ou negatrf apres transformation par l'op6rateur moment cin6tique selon z. Dans ('expres- 
sion I on a pos6 : 

a 

j(Wx + jw y ) = w 1 e!0 

15 

ou Wi est ('amplitude de I'excitation et ou 0 est sa phase d'emission par rapport aux axes x ety. Cette phase 
est conventionellement nulle quand remission est selon I'axe - y. 

Pour integrer liquation I par rapport au temps on procede d'une maniere pr6f6ree a une approximation, 

20 de type connu, dite par "impulsions dures" ou encore "hard pulse". Selon cette approximation on considers 
que les protons sont soumis, pendant I'excitation, a deux actions differentes. Une premiere action permanente 
est relative a la presence du gradient de champ orientateur dans la tranche 7. Elie a pour effet le decalage de 
pulsation de resonance w. La deuxieme action Wi correspond a I'excitation proprement dite. Dans la pratique 
Wi evolue dans le temps et c'est son evolution que Ton recherche. En echantillonnant on peut la representor 

25 sous la forme d'une suite de valeurs bloquees en paliers pendant des p6riodes eiementaires et evoluant bru- 
talement, d'un palier a un autre, de la fin d'une p6riode au debut d'une autre, (.'approximation evoquee consiste 
a remplacer Involution par palier de la deuxieme action par une succession d'impulsions de type equivalant a 
une fonctlon de Dirac, de duree s tres courte et dont I'amplitude W| est egale a w t (t) - x/e ou Wi (t) est la valeur 
des paliers de I'excitation radlofrequence et ou x est la dur6e d'une p6riode 6l6mentaire. L'effet de cette approxi- 

30 mation est de simplifier notablement liquation I. En effet on peut dire que pendant les impulsions dures, elie 
s'6crit : 




parce que W| est tres superieure a l'effet du gradient de champ orientateur. 

Hors des p6riodes d'appiication des excitations, cette expression peut s'ecrire simplement 



45 
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dt 
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dt 
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55 puisqu dans ce cas I'impulsion d Dirac n'est pas presente. On peut integrer chacune des deux expressions 
pr6c6d ntes par rapport au temps. L calcul stclassiqu , et on peut ecrire respectivement : 
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25 



Dans ces expressions les amplitudes de probability sont affectees du signe - ou du signe + selon qu'elles 
correspondent a ces probabilites avant ou apres Implication de I'impulsion de Dirac i correspondant. Eiles sont 
affectees d'indices i ou I + 1 selon le rang i de I'impulsion de Dirac qui les precede. Dans ces expressions 0, = 
W| - s = Wi (t) • x. Si par ailleurs on pose pour simplifier les expressions : 



e>* = z 
cos e/2 = c, 



ou 



et 



4 



J 

e sin e./2 = s. 



35 



on peut ecrire qu'apres une n ieme impulsions de Dirac, en omettant les signes 



40 



n 



n 



c - s 
n n 



s c 
n n 



z 1 ' 2 0 



0 2 



-1/2 



n-1 



'n-1 



CII) 



45 



De cette ecriture on va deduire des expressions de a n et b n en fonction de a^ et b^. Et dans ces der- 
nieres expressions, en remplacant a n et b n par : 
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et 



on peut ecnre : 



a = 



n/2 



V 2 



n/2 



(III) 



10 



15 



— -1 - 

a = c a ,-S2 6 
n nn-1 n "n-1 



s_ a . + c z " 1 6 
nn-1 n Tn- 



1 



(IV) 



Si on suppose connue 1'orientation du moment magnetique avant u ne premiere impulsion de Dirac, en par- 
tlculier s'il est orients parallelement a B 0 , on peut dire que : 
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et done 
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Dans ia demfere expression, a<, et p 0 dependent de Cq et So qui sont repr6sentatifs de ('amplitude de la 
premiere impulsion de Dirac numerates 0. Cette formulation peut encore s'6crire : 



40 



et 



45 



0-0 

Po = s o = So 2 



so ou z° vaut conventionnellement 1. Apres la deuxieme impulsion (numerates n = 1) on peut Scrire 



55 



a i 8 c i c o s o z ~ l - h o 2 ~° + h J 2-1 
*i = s i c o ■ c i s o 2-1 - 4 2 "° + s{ 



De meme apres une troisieme impulsion (numerate n = 2) on peut ecrire : 
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°2 * C 2 (C 1 C 0 * S l S 0 2 " l) * V~ 1(S 1 C 0 * C l S 0 2 " 1} 
soit a 2 = h| z" 0 + h 2 . z~ 1 + z ' 2 . 

et Ton pourrait 6crire de m6me pour : 

a -„ 2 ,-° 2 -1 2-2 

% 

ft 

Alnsi de suite et plus g6n6ralement aprds ia n plus unifcme impulsion, num6rot6e n, an et p n s'6crivent 

a n " h 0 2 + h i z + - + h n 2 
Pn ~ g 0 2 + «l 2 + - + S n 2 
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Avec cette expression on peut faire deux remarques importantes. Premidrement a et p reprise ntent, k un 
changement de variable prfes, les amplitudes de probability de la fonction d'onde transform^ par I'opSrateur 
moment cin6tique selon I'axe z. On peut d'ailleurs montrer, qu'au mfime titre que le carr6 du module de a addi- 

25 tionn6 au carr6 du module de b vaut 1, que le carr6 du module de additionn6 au carrt du module de p n vaut 
1 6galement. En effet la probability pour que le moment magndtique existe vaut toujours 1. Deuxfemement on 
remarque que ('expression de a n et p n est une expression polynOmiale analytique en z. Or z correspond k la 
pulsation de resonance en fonction du lieu (k cause du gradient) des protons concem6s dans I'espace excite. 
Done z est dans I'expression du moment mag n6tique, le terme repr6sentatif de la frequence de resonance. Au 

30 meme titre que cq ei sq eiaient deienriines par "amplitude de ia premiere excitation on peut dire que !es coef- 
ficients h? et ffi $ont d6termin6s par la succession des differentes impulsions, jusqu'S cede de num£ro n. Ces 
coefficients sont cependant exprim6s, eux, indSpendamment du lieu des protons concern^ e'est k dire ind6- 
pendamment de ta frequence. Autrement dit, la connaissance de tous les hj et gy permet de d6duire failure du 

moment magn§tique k la fin de ('excitation quelle que soit la frequence de resonance, e'est k dire quelle que 
soit la place des protons dans la tranche. 

Dans ('invention on se prtoccupe alors de rechercher de tels coefficients et g n { en imposant des contrain- 

tes li6s aux valeurs de w. En quelque sorte, on va imposer que les expressions de otn et p n soient les Equations 
d'un filtre ayant une vateur donn6e dans sa bande passante et ayant une autre vateur donnSe k I'exterieur. 

On peut 6crire en reportant la valeur an - 1 et de p„ _ ^ en fonction des coefficients h n - 1 et g n 7 1 de I'expres- 
sion V dans I'expression IV de a„ et p n , eux-mSmes Merits en fonction des coefficients h? et g; : 



. n . n - 1 - n - 1 

n. = c • n. - s • e. , 
i n 1 n 6 l - 1 



et n .n - 1 n - 1 

e. = 5 . n. + c • g. . 
*i n i n °i - 1 



CVI) 



50 L'expression IV peut en outre fitre invers6e et devenir : 



a . = c a ♦ s ft 
n-1 n n n r n 



B . = z (- s a •+ c 6 ) 



CVIIIJ 
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10 



15 



A ce stade du calcul mathgmatique on peut faire les considerations suivantes. Si on s'est donn6 a,, et p n> 
repr6sentatifs d'un moment magnetique a la fin d I'excitation, on peut, connatssant Cn et s„ r presentatifs de 
la dernier impulsion de Dirac de I'excitation, dSduire les valeurs du moment magnetique, a,,- 1 et 0 n _ t, juste 
avant Implication de cette impulsion. De proche en proche, en determinant a chaque fois Fimpulsion prece- 
dents on peut remonter revolution du moment et determiner, a I'envers, ITiistorique de I'excitation depuis son 
origine jusqu'a sa ftn. La question est alors de savoir si on peut choisir c n ou s n comme on le veut ou au contraire 
si ces valeurs sont imposees. Que ces valeurs soient imposees revient a dire que impulsion numero n est 
connue en amplitude. En examinant VIII et en se rappellant que a n et p n sont des polynfimes de degr6 n en 
z- 1 il vient Imm6diatement : 

0 = c n + s £ 
n n n e n 

et 0 = - s nV c n 8 0 



20 



En effet le polyndme de degr6 n - 1 en r 1 repr£sentatif de a n _i ne peut contenlr des termes de degr6 n 
en z- 1 . De meme le polynome de degre n-1 en r 1 representatif de p n _i ne peut comporter des termes en 
z~ 1 . Ce couple de conditions revient a poser les questions suivantes : 



25 



s e' ^n. sin © L 

n n 12 

ST cos e 12 

n n 



t 8 

2 = 



n 



g 



n 



S 



n 
0 

n 



n 



30 



Ceci revient a imposer que © n et 0 n soient les deux inconnues d'un systems de 4 equations a deux in con- 
nues. En principe on devrait se heurter ici a une impossibility de resolution. 

Mais le fait de dire que la probability que le moment magnetique existe vaut 1 revient a ecrire que : 



35 



so 



(VII) 



ou encore : 



a n a n + p„p„ = 1 (VII) 

40 Cette expression VII peut etre developpge en remplacant a n et p n par leur valeur don nee dans F expression 
V. En effectuant ce dernier calcul et en classant les coefficients obtenus selon les puissances de z auxquelles 
ils s'appliquent on peut 6crire : 



45 



55 
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h 0 h n-l + h l h n + «0 «n-l + «l *n 



s 0 termes en z 



-n 



= 0 termes en z 



-<n- 1) 



•••• 



•••• 



= 0 



10 



h 0 h 0 + h l h l * 



,+ h n + 
n n 



15 



n n n n 
8q S 0 + g l g l + 



n n 



a 1 termes en z 



-0 



= 0 



20 



h 0 h n 



n 



+ Sq 



S 



n 



= 0 termes en z 



25 



30 



35 



40 



45 



50 



Ce systeme d'gquations comporte n + 1 Equations, et non pas 2n Equations, & 2n + 2 Inconnues. En effet 
les Equations dont !es termes sont en tt* ou en z 1 sont les mimes. Par contre i! y a bien 2n + 2 termes en hy et 

ffl qui sont inconnus. On remarque que la premiere et la derniere equations reduisent l'impossibilit6 de deter- 
mination de 0 n etde 0 n telle qu'elle apparaissait precgdemment. En effet les deux derniers termes de Pexpres- 
sion IX sont identiques entre eux selon la condition imposee a la probability de I'amplitude du moment 
magnetique. En consequence les deux Equations a. deux inconnus © n et 0 n sont solubles. II en r6sulte que 
0. et 0* sont fixes par ia nature des polyndmes en hj et en g? qu'on a retenus. On peut done calculer a n . 1 et 

p n _ 1 en application de I'expression VIII. D'ou on deduit 8 n - 1 et 0 n _ 1a Ainsi de suite on peut calculer la sequence 
des impulsions qui composent I'excitation. 

A ce stade de I'explication on peut risumer la situation en disant qu'on peut choisir deux polyndmes de 
degri n parmis un ensemble important de polyndmes possibles (a. cause du systeme de n + 1 equations k 2n 
+ 2 inconnues). On sait, une fois que de teis polyndmes ont ete choisis, qu'on est capable de calculer la 
sequence des excitations sous formes d'impulsion de Dirac qui ont conduit, d paitird'un moment magnetique 
d'orientation donnee, aux polyndmes en question. En refaisant aiors {'approximation par impulsion dure, mais 
& i'envers cette fois, on peut donner une vaieur pseudo continue ou mime continue a i'ampiitude de ('excitation. 
Autrement dit on vient de transformer {'expression physique du phenomena. Sur cette derniere expression on 
peut effectuer simplement des optimisations. Toiites les optimisations porteront sur un choix, en fonction de 
criteres d'optimisation donnes, des coefficients h? et gj dont les vafeurs sont conformes aux conditions enon- 

cees dans I'expression VII. Les calculs d'optimisation en outre reviennent & des calculs de filtres ftequentiels. 
On sait imposer des conditions d'amplitudes et de phases & de tels filtres. 

D'un point de vue macroscopique I'effet global de I'excitation radiofrequence combined au mouvement de 
precession est de soumettre I'aimantation a une rotation de parametres angulaires ©(w), psi (w) et phi (w), ou 
B psi et phi sont les angles d'Euler definissant cette rotation, w est la pulsation de resonance consideree. Cha- 
cun des trois angles Q ps! et phi est fonction, comme a et b, de I'historique de I'excitation W1 (t) et peut dtre 
trouvi, pour chaque w, par integration sur le temps des equations de Bloch. II existe une conrespondance 
directe entre les angles d'Euler de cette rotation et les amplitudes de probability, a et b, solutions de I'equation 
de Schrodinger (Albert Messiah, "Mecanique quantique", Tome II, pages 918, 919 et 922). Ces corresponden- 
ces s'ecrivent : 



a = £ exp (j phi/2) cos (9/2) exp (j psi/2) 



55 



b = e exp ( - j phi/2) sin (9) exp (j psi/2) 



10 



• t 



EP 0 316 312 B1 



10 



25 



on 



avec £ = * 1 



On peut inverser ces relations et obtenir : 

cos 0 = (aa- bb) = (eta- pp) 

soit 

cos 9 = 1 - 2 pp = 2 aa - 1 

et 

15 psi = Arg(a • b) = Arg(z n • a • P) 

phi = - Arg(a . b) = - Arg(a • P) 

On peut chercher a r6soudre le problems consistant a approximer, sur tout un domaine de w, un angle de 
20 nutation donn6 0 (w) sans s'interesser aux angles psi et phi. 

On salt que si ('excitation est un train de n impulsions de Dirac alors ies a etg resultants sont des polynOmes 
d'ordre nenr 1 . II est facile de voir alors, compte tenu du fait que |z| = 1 et done que le conjugue de z est £gal 
a z -1 , que pp et aa sont des"polyn6mes" de "degr6" 2 n (puissance de z de - n a + n). Ces poIynGmes ont par 
ailleurs la particularity d'etre sym6triques et de pouvoir s f 6crire : 



i 



n - i 

Trouver Ies coefficients p r qui approximent au mieux une reponse en frequence (cos 0 (w) - 1)/2 donnee 
35 est un probleme dassique de filtrage. Voir par exemple A.V. Oppenheim, R.W. Schaefer, Digital Signal Pro- 
cessing, Chapitre 5. On utilisera de fagon preferee Talgorithme dit d'6change de Remez expose dans cette refe- 
rence et publlfe dans Ies IEEE Transactions on Circuit Theory, CT-19, 189-194, Mars 1972 par MM T.W. Parks 
et J.H. McClellan sous le titre "Cheby s hev approximation for non recursive digital filters with linear phase". 
En plus de son caractere d'optimalite cet algorithme presente Tavantage de pouvoir rendre compte facilement 
40 de la contrainte : 

0<pp<1 

Une fois Ies p, determines on peut trouver Ies polynfimes p(z). En effet si z, est un z6ro de p alors p(z) = 
45 pp a comme zero z, et 1/Z|. A Tinverse connaissant les 2n _z£ros de p(z) on construit p en choisissant parmis 
chacun des n couples (z,, 1^,) un de ses elements z, ou 1/Z|. II y a done n choix binaires possibles et, pour un 
PP donne, 2 n solutions possibles pour II en va de meme pour a : pour aa donne il y a 2 n solutions a possibles. 
II y a done, pour un cos 9 (w) approxime au mieux par un p(z), 4 n solutions pour les 2 autres angles d'Euler 
psi et phi. 

so On peut aussi, par recherche des zeros de p(z) et 1 - p(z), examiner toutes les solutions psi et phi, mais 
on conpoit que cette operation est longue et fastidieuse . En fait le nombre de solutions possibles n'est pas aussi 
grand qu'il y parait de prime abord car Falgorithme de Remez tend a placer le plus de zeros possibles sur le 
cercle unite |z| = 1 , et dans ce cas z, = 1 H v 

En ce qui conceme a on peut general m ntserestreindre a chercher une solution a phas minimum. Cette 

55 solution correspond a ch tsir entre ies zeros Zj 1 1/|z, celui qui est a I'interieur du cere) unite. En fait cette 
solution peut etretrouve sans avoir a chercher les z6 no de 1 - p(z) grfice a un alg rithme simple inv ntepar 
MM. R. Boite et H. leich, publie dans Signal Processing, 1981, pages 101 a 108 par N rth Holland Publishing 
Compagny dans un article d nt le titre est : 'A NEW PROCEDURE FOR THE DESIGN OF HIGH ORDER MINI- 
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MUM PHASE FIR DIGITAL OR CCD FILTERS". D'un point de vue physique la soluti n a phase minimum pour 
asejustifie pari faitqu'ell c nduit, apres inv rsion des equations d Shrodmger par IX, a I'energie de I'exci- 
tation la plus faible. En eff t nsaitqu pannis toutesles solutions a de aa= 1 -p, la solution a phase minimum 
est celle qui a le coefficient d'ordre z6ro h n 0 le plus grand. Or en reprenant les relations VI on trouve que : 

A 

. n cos ' n-l # cos cos 0 

h 0 s 

10 Si Ton prend le logarithme de cette valeur, et si on considere que les contributions apportees par chacun 
des W| a la nutation totale du moment magnetique sont faibles, on peut ecrire, apres avoir effectue un deve- 
loppement limite de cosinus 9, : 



15 



50 




n 



Log hj - -|X<v : 



Comme on a vu que hj etait le plus grand on en d6duit que la somme des carr6s des 8j est la plus faible. 

2° || en resulte que les impulsions, dont la succession provoque une somme des deux angles d'Euler psi + phi 
minimum, sont justement celles qui correspondent a I'energie minimale d'excitation pour un angle de nutation 
0 donne. Ce resultat n'etait pas connu precedemment 

Une approche plus restrictive de la definition des polynomes a et g consiste a imposer des le depart des 
contraintes sur les phases arg (a) ou arg (b). En particulier on peut imaginer que arg (b) = 0. Ceci est equivalent 

25 a imposer que les coefficients du polynfime £ soient symetriques, c'est a dire que 

n n 
8 55 & i 

30 n - i 

Le polynflme ^ peut alors 6tre trouve par optimisation de Tangle de nutation 9 (w) grace a Talogrithme de 
Remez. On en deduit alors pp et aa. On choisit alors a a phase minimum. La figure 2 montre i'allure d'une 
impulsion a 180° trouvee en respectant ces conditions. Cette figure est associee a un premier tableau, TAB 

35 1 1 de valeurs de Tamplitude de I'excitation en fonction d'un certain nombre de periodes (les periodes sont don- 
nees en temps reduit et les amplitudes en radian s). Parce que la forme de cette excitation est symetrique dans 
le temps, seules les 1 + n72 premieres periodes sont tabulees. Dans Texemple H y en a 59. Pour des raisons 
technologiques cette impulsion peut conduire a des puissances cretes, au centre, trap fortes pour les genera- 
teurs 3 utilises. On peut alors modifier son allure en prenant une des autres solutions de a. On aboutit alors a 

40 une impulsion dont I'enveloppe est celle representee sur la figure 3, et qui correspond, dans les memes condi- 
tions, a un deuxieme tableau de valeurs : TAB 2. A Texperience il apparait que le respect a 3 ou 5% pres des 
valeurs indiquees dans ces tableaux n'est pas prejudlciable au resultat des experimentations de RMN qui les 
utilisent 

45 



55 
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TAB 1 



TAB 



-3. 9049108 88 20021E-002 
- 6. 9663 1 9801 34976E-003 
-9. 23 1 66S98398407E-003 
-9. j 967746&O407O8E-0O3 

- 8 . 8 34 52823996 1 23E-003 
-6. 12391456332324E-O03 
- 7 . OS94 1 632094723E-0O3 
-5. 6347 7337 2383O7E-O03 

8623394 10 13924 E-003 
- 1 . 7e46*4T,2l384SBE-O03 
6. 25474B38024969E-OO4 
Z ■ . Z 64 40690*4368 ^»003 
6 . 0 93830 7459 1 332E-OC3 
9 . 032 16904 39039 3E-O03 
1 .203351 42S41780C -002 
1 .50 19*7681 9187 4E-002 
1 . 783 1 74 1 35947S7E-002 
2 . 0433734 637 4239C -002 
2 . 27244 1 6828332 1 E-002 
Z'. 436301 46934734E-002 
2 . 387 1 028 7 1 53642E -002 
2 . 4357 6B6 1205690E-O02 
r.6543135S434137E-002 
2. 3759] 1 39192740E-002 
2.41343403133147C-002 

2. J 69349 1 7603766E-O02 

1 . 8371 3380632749C-0O2 
1.41 9430421 093O3C-O02 
9. 264S7660883S27£^003 

3. 46763360448094C-O03 

- 2. 92329 ] 0073399X-003 
-9. 84733906633224C-O03 
- 1 . 7 1 S30 1 737844S8C-002 
-2. 4637133381 6464E-O02 
-3. 2 144221 94442 14C-O02 
-3. 9349381 40166O2E-OO2 

•-4 . «042 1 240029 1 87C -002 
-3. 19305174098276C-OO2 
-3. 6636469474841 0€ -002 
-5. 98922285223749E-O02 
-6. 13121 31 9339 390C-002 
-4. 03981 140770231C-O02 
-3. 744947391 12293E-O02 
•5. 1 3*20362699079* ~O02 
-4. 26429999744443C-O02 
-3. 0305368 1 130646E-O02 
-1 • 48418 1085734O1E-OO2 

4. 48331 838317662C-003 

2. 773 14 167977363C-O02 
3.3G441462341713E-002 

0. 4497290324*637C-002 

1 . 2193136I422393C-001 

1 . 40897043098379C-O01 
2. 0220306397 1 342E-O01 

2. 43859077421 071C-O01 

2. 828239 39221 93 lE-OOl 

3. 151 8*97921 8S70C-001 
3. 368287B0039443E-OO1 
3.44S191O1220646E-O01 



-4 . 4 34 20 1 332 1 999 1 £-002 
-2. 22S06739796626E-002 
*2.4984902O943l08E-O02 
-2. 6899O157438880E-002 
-2. 77707348 036423C-OO2 
-2.7359471 191431 1 E-002 
-2.34479434337933E-002 
-2. 1 863347 24 M 27 9E -002 
-1 . 44933443410009E-002 
*9 * 3973 1 1 644 45 1 4 1 E-O03 
-3. 1 9 1 32207744 244E-0O4 
1 .030290774 08931 E-O02 
2.23249900304281E-002 
3 • 339084272 37378E -002 
4 . 90274 190413313C-00? 
4.26431 148417848E-002 
7.33446371967290C-002 
• . 494792281 31 243E-002 
9.4015522I41S1 16E-002 
1 « 01921 802249669E-00I 
1 • 039 1 3 1 38439 1 47E-001 
.1.013801 87889778C-O01 
9.391 3332 1 324 1 70E-002 

0. 1 3443343002227E-O02 
4. 3739022 14774O7T-O02 
4 . S 379O403273033C-O02 

1. 474845481 49141E-002 
- 1 . 948392344001 07E-O02 
-4. 844244 1014884 1C-O02 

2?ftS2S3l 4B6asX-ttl2 

-1.1 7423798321 197C-O01 
-1.30471 443440478C-OO1 
•I . 79244867442907E-001 
-2. 0143901384 24938 -OOl 
-2. 1 47492834 13040C-001 
-2. 1 743478448 1U7E-001 
-2. 067 149 1 89O0338C-O01 
-1 . 0903O077064343E-O01 
- 1 . 99738970370414E-001 
- 1 . 22847407684439E-00 1 
-8. 040574 1 2979931E-002 
-3. 43898608033487C-002 
1 . 34844833844243C-O02 
4* 14742832514041C-O02 
t . 08382404007442E-001 
1 .3237236040 1,07 1E-001 

1 . 908948709830O3C-OO1 
2.21041 11222O419E-001 

2. 431 14747839940C-O01 
2. 3339O344178343C-001 
2.32374007632722C-001 
2.41 98 3370764 297E-001 
2. 2408370931 1 999C-001 
2. 01994161 173123E-001 ■ 
1 • 78903 147399749E-001 

1 • 57431 97649*98 IC-OOI 
1 . 40483393387081E-001 
1 . 29003722671 104C-001 
1 .24 1 64 37623286 1E-001 
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Revendicatl ns 

1. Procede d'excitation par un signal de radiofrSquenc , dans un appareil d'experimentation a resonance 
magnetique nucleaire comportant des moyens electroniques de traitement de I'information, excitation en 

5 reponse a laquelle (es bascuiements de moments magnetiques de spins nucleates contenus dans un volume 
d'investigation determine soumis a ('experimentation sont definis en amplitude et en phase pour un ensemble 
de frequences de resonance donne, ledit procede etant caracterise par ies etapes suivantes : 

— on introduit dans lesdits moyens ('equation de Schrodinger relative a I'amplitude de probabilite de la 
fonction d'onde de spins de ces particules resonnant a une frequence de resonance donnee, rhamiltonien 

10 deladite equation etant constitue par une fonction dependant de la frequence etde la forme d'une impulsion 

d'excitation appliquee a ces spins, 

— on integre par rapport au temps ladite equation de Schrodinger dans lesdits moyens et I'on en deduit 
Ies expressions analytiques polyndmiales des amplitudes de probabilite a la fin de ('impulsion en fonction 
de la frequence de resonance, ces expressions etant d'une part, de degre egal a un nombre de periodes 

15 elementaires dont la somme des durees est 6gale a la duree de ladite impulsion et, d'autre part, affectees 
de parametres constitues par des coefficients representatife du profit d'excitation effectif cause par ladite 
impulsion et subi, pendant lesdites periodes, par Ies spins aux differents points dudit volume d'investiga- 
tion, 

— on optimise par lesdits moyens lesdits coefficients de facon que lesdits polynomes s'approchent aux 
20 mieux d'une forme pr6detenminee d6siree, 

— on en d6duit, pour chacune desdites p6riodes 6l6mentaires, Ies caract6ristiques d'un signal d'excitation 
conduisant a ladite forme pr£d6terminee desir6e, 

— on excite Ies spins dudit volume avec ce signal d'excitation. 

2. Proced6 selon la revendication 1, caracteris6 en ce que pour integrer on effectue une approximation 
25 par impulsions dures. 

3. Proced6 selon Tune quelconque des revendications 1 ou 2, caracteris£ en ce que pour optimiser on 
recherche Ies expressions polyndmiales a dephasage minimal qui correspondent a une excitation la moins 
puissante possible. 

4. Proced6 selon la revendication 3 f caracteris6 en ce que pour exciter on excite avec une excitation dont 
30 Failure est celle de la figure 2 representative, du tableau suivant, ies periodes eiant donnees en temps reduit 

et le amplitudes en radians : 
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TAB 1 



8 . <>6*3l 9801 34 976G-003 

9. r3J6659eS98407E-003 
19677465040708E-CO3 

e. 834 32825998 123E-003 

6. 12591456552324E-003 
7 . C596 1 652094723E-OO3 
5.63477557258307E-003 
3. 8625391 1015924E-003 
1 .7c46*4r.Sl3B43BE-003 
6. 25476838024989E-0O4 
2*440690345683.6-003 
6 . 095630 7659 1 332E-OC3 
9 . 0321 6904390593E-O03 
1 , 20555 1 4236 1 780E-002 
1 .5019A768191874E-002 
1 .785174 1 3396737E-002 
2. 0453734637 *23?E-002 
2. 27246 1 68283321E-002 
S« 456301 46956734E-002 
2. 587 1 028 7 1 3364 2E -002 
I'. t»357 G66 1 203690E-002 
2.6S4313S5636157E-O02 

1. 3759 1 1 39 1 92740E-002 
2.41 345403 1531 4 7E-O02 

2. 1 69569 17803766E-002 
1 . 837 1 3380832749E-002 
1.41 9430621 09303E-O02 
9. 2045766O88S527E<*OO3 
3= 467^5560648096E-003 

- 2 . 92329 1 00733593E-003 
-9. 8473390B633224E-OO3 

- 1 . 7 1530 1 73784653E-002 
-2. 4637133381 66^4E-O02 
-3. 214422194*4214E-O02 
- 3 . 936938 1 401 6602E-002 
•-4 . o0621 260029 1 37E-O02 
•5* 19303174898276E-O02 
-3.66364696748610C-002 
-5. 98922285225749E-002 
-6. 13121319539390E-002 
• 6 . 0598 1 1 4077 823 1 E-002 
-5.744947391 12293E-002 
-5. 13620362699079E-OO2 
-4.26629999764663E-002 
-3.05053681 138646E-002 
-K48618lOe573401E-002 
4. 48331 8383l7^62E-003 
2. 773 1 4 1 67977343E-002 
3. 504 6 1 4 6234 17 1 5E-O02 
8. 64972905264657E-002 
l«2193156l422393e-001 
I • 60697063098379E-O0 1 

2. 02203063971342E-OO1 
2.43859077421O71G-O01 
2.82823939221931E-O01 

3. lSie69792ie870E-001 
3. 36828780839443E-O01 
3.4451910I220846E-O01 
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5. Proc6de selon la rev ndicatfon 3, caracteris6 en ce que pour exciter on excite avec une excitation dont 
Taliure est c lie de la figure 3 representative du tableau sutvant, les period s etant donnees en temps r6duit 
et les amplitudes en radians : 
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TAB 2 



•6. 436201 532 1 999 I E-002 
-2« 22S06739796626E-002 
-2.4986902O943108E-002 

-Sr. 6e99oit743eesoe-oo2 

-2.77707368058623E-0O2 
-2. 7359471 19! 651 1 E-002 
-2,544 794 34537933E-002 
-2. J663S472488279E-002 
-1 . 6493346S61OOO9E-002 
-9.29751 16646S141E-003 
-3. 19132207764246E-OO4 
1 .030290774 08931 E-002 
2. 2324950030628 1 E-002 
3 • 53906627237376E-002 
« . 90274 1 906 1 53 1 5E-002 
6. 2643) 1 484 17868E-O02 

7.55466371967290E-002 
8 . 69479228 151 243E-002 

9.601 552214 15U6E-002 
1 •01921802249669E-001 
1 . 039 1 3 1 384 39 1 67 E-OO 1 
.1 .01 3801 87 88977 8C-O01 
9. 39 1 3532 1 524 1 78E-002 
6. 13463563O02227E-O02 
6 . 37 39022 1 6774 07E-002 
4 . 1 37 90W3273033E-C02 
I. 474865481 491 61 E-002 
- 1 . 54 8392344001 07E-O02 

-4 . 644244 1 01 6884 1 E-002 
-0. 2985233 1 664 B32E-002 
-1.1 762579552 1 1 97E-00 1 
•1.504714636486786-001 
-1 .79244 B67442907E-00! 
-2. 0 1 43901 3642493E-O01 
-2, 147492e5613060€-O01 
-2. 17434764481 1 17E-001 
-2 . 067 1 69 1 8900358E-001 
-1 • 690340? 7066345E-001 
-1 .3973897037061 6E-001 
- 1 . 22&47407666439E~00I 
-8 • 04 057 6 1 297993 1 E-002 
-3. 4389&606033687E-0O2 
1 . 34 844 63364 624SC-002 
6. 167 4 26525 14041 E-002 
1 . 063826O6007642E-001 
1.3237256O401O71E-O01 

1 . 9O696870965O03E-001 

2. 21 84 1 1 1222041 9E-001 
2.431 14767633960E-OO1 
2.33390364178345E-001 
2. 525766O7632722E-0O1 
2.41 983570764 297E-001 
2. 24 0837053 1 1 999E-001 
5.01994161 175123E-001 
1.78905167599749E-O01 
1 . 37651 976499981E-0O1 
1 .4O-*853933870*8lC-0Ol 
1 . 2900372267 1 1 06E-O0 1 
1 . 24 1 64376232861E-O01 
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Pat ntanspruche 

1. Verfahren zur Anregung durch ein Hochfrequenzsignal in einem NMR-Untersuchungsapparat mit elek- 
tronischen Infomnationsverarbeitungsmitteln, wobei die durch die Anregung hervorgerufenen Umkiappungen 
der magnetischen Momente der Kernspins, die in einem der Untersuchung unterworfenen, vorgegebenen 
Untersuchungsvolumen enthalten sind, fur eine gegebene Gruppe von Resonanzfrequenzen in bezug auf die 
Amplitude und die Phase bestimmt werden, wobei das Verfahren gekennzeichnet ist durch die folgenden 
Schritte ; 

— Eingeben der Schrodingergleichung, die auf die Wahrscheinlichkeitsamplitude der Spin-Wellenfunktion 
dieser bei einer gegebenen Resonanzfreqeuenz Resonanz zeigenden Teilchen bezogen ist, in die Mittel, 
wobei der Hamilton- Operator der Gleichung durch eine von der Frequenz und der Form eines an diese 
Spins angelegten Anregungsimpulses abhangige Funktion gebildet wird, 

— Integration der Schrodingergleichung bezQglich der Zeit in den Mitteln und Abieitung der anaiytischen 
Pofynomausdrucke der Wahrsheinlichkeitsamplituden am Ende des Impulses in Abhangigkeit von der 
Resonanzfrequenz, wobei diese Ausdrucke einerseits von einem Grad sind, der gieich einer Anzahl von 
Elementarperioden ist, deren Gesamtdauer gieich der Dauer des Impulses ist und andererseits mit Para- 
metern versehen sind, die durch Koefflzienten gegeben sind, die das effektive Anregungsprofii darstellen, 
das durch den Impuls bewirkt wird und wahrend der Perioden mitteis der Spins an verschiedenen Punkten 
des Untersuchungsvolumens anliegt, 

— Optimierung der Koeffizienten durch die Mittel, derart, daR sich die Poiynome am besten einer 
gewunschten, vorgegebenen Form annahem, 

— Abieitung der Eigenschaften eines zur gewunschten, vorgegebenen Form fuhrenden Anregungssignals 
aus den optimierten Koeffizienten fur jede der Elementarperioden, 
— Anregung der Spins dieses Volumens mit diesem Anregungssignal. 

2. Verfahren gemaB Anspruch 1, dadurch gekennzeichnet, dad zur Integration eine Approximation mitteis 
kurzer Impulse ausgefuhrt wird. 

3. Verfahren gem5B einem der Anspruche 1 oder 2, dadurch gekennzeichnet, da& zur Optimierung die 
Polynomausdrucke mit minimaler Phasenverschiebung gesucht warden, die einer Anregung mit der gerings- 
tmoglichen Leistung entsprechen. 

4. Verfahren gemaU Anspruch 3, dadurch gekennzeichnet, daR zum Anregen mit einer Anregung angeregt 
wird, deren Veriauf derjenige der Fig. 2 ist, die die folgende Tabelle darstellt, wobei die Perioden mit verkurzter 
Zeit und die Amplituden in Radian gegeben sind : 
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Tabelle 1 



3, 90491O8882O621E-OO2 

8 . *60M 9601549 76E -003 

9. 2Z J 66396298407E-OO3 
1 9*? 7 4 6304 07 08 E -003 

8. 834}2823999l23E-003 
6. 12391436332324E-003 
7 . 0396 1 6520947236-003 
2. 63477357236307E-003 
3.6622391 1013924E-0O3 
1 .7c4G6432l38436E-C03 
*. 22476638024 989E-004 
«6440690345683£-003 
6 . 0938307639 1 332E~0C3 
9 • 032! 690439O393E-OO3 
1 • 20323 1 4236 1 780C-002 

1 . 30 1 9 a7 68X92 874 E~0O2 
1 .7831741 3396737E-002 

2. 0433734637 6239E-002 
2. 272461 68283321E-002 
2. 436301 46936734E-002 
2* 387 1 028 7 1 33642E-002 

63376S61203690E-002 
2. 6343 1 333636 1 37E-002 
2. 57 29 U 39 1 9274OE-0O2 
2* 4 1 343403 1 33 1 47E-002 

2. 1 69369 17603766E-O02 
t . 837 1 33806 J2749E-002 
1 • 41 943062 109303E-O02 
9.2043766O683327£<*OO3 

3, 46763366648096E-OO3 
■2.923291 0O733S93£-AA3 

•9i 84733906633224E-OO3 
•1.713301 73784653E-002 
•2. 46371333© 1 M44E-002 

3.21442219444214E-O02 
>3.936938t 401 6602E-0O2 
•4 . o062 1 260029 1 37E-0O2 
3. 19302174898276E-O02 
-3. 6636 4 696 74 86 1 0E-O02 
•3. 98922283223749E-002 
•6. 131213 19339390E-O02 
6.03981 14077 623 1E-O02 

•3.744947391 12293E-002 
•5. 1 36203 62699079E-O02 

26629979764663C-O02 
•3. 03033681 138644E-002 
• 1 . 486 1 8 1 0637340 1 E-002 
«. 48331 838317662E-O03 
2* 773 141 679 77563E-O02 
.3. 3046 1 4 6234 171 3C-O02 

G.64972903264637E-002 
U2l9313«1422393€-00l 
U 60e97063O96379E-O0l 
2. 0220306397 I 342E-O01 
2. 438390 7742 1071 E-OOl 
2* 82823939221931E-O01 
3. 12166979216670E-OO1 
3. 368267eo639443E-OOl 
3.442191 1220846E-OQX 
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5. Verfahren gemafc Anspruch 3, dadurch gekennz ichnet, dad zum Anregen mit einer Anr gung angeregt 
wird, deren Verlauf der jenige von Fig. 3 ist, die die folgende Tabeiie darstellt, wobei die Perioden mit verkurzter 
Zeit und di Amplituden in Radian gegeben sind : 
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Tabelle 2 



-6 . 4 Zb 20 \ 332 1 999 1 £-003 
-2. 22S06739796626E-O02 
-2.4986902O9431O8E-O02 
-2. 68990157438880E-OO2 
-2 . 7 7 7 07 36 8038623C -002 

-2.7339471 191631 1E-002 
-2.544794^537?33£-oo2 
-2. 1 663547248827 9C-O02 
-1 . 649354656 J 0OO9E-O02 
-9. 29 7 5 1 1 664651 4 1 E-O03 
-3. 1 91 32207764246E-OO4 

1 . 03029 0774 0693 1 E-O02 

2. 23249500506281 E-002 
3 • 539 08627 237378E-002 
* . 90274 1 906 1 55 1 SE-002 
6.26451148417868E-002 
7 .554 6637 1 96729 0C-O02 
6* 69479228 151 243E-002 
9. 601 552214 15116E-002 
1 .01921 802249669E-OO1 
1 • 039 1 3 1 384 39 1 67E-00 1 
1.01380187889778E-O01 
V . 39 1 5552 1 524 1 7 6E-O02 
©. 134 63565002227 E-O02 
6.37390221677407E-002 
«. 1379060327 3033E-O02 
1 • 4 74 8654 81 49 1 61 E-002 

-1 .548392344001 07 E-002 
~4 . 644 244 1 016884 1 E-002 
-6 . 29 65233 1684 652£~002 
-1 . 1 7*25793321 117E-O01 
-1 • 30471465648676E-001 
-1 . 79244867442907E-001 
-2. 01 439013842493C-OO1 
-2.1 47692636 13O60C-O01 
-2. 1 7434784481 1 17E-001 
•2 . 067 1 69 1 8900358E-00 1 
•1 • 8903007 70*634 3E-00 1 

-1.39738970370616E-001 
- 1 ♦ 22847 407686439E-00 1 
-6.O4O57612979931E-O02 
-3. 43898608033687E-002 
1 . 34 8448338 4 6245E-002 
6. 16742852514041E-002 

1 • 06S82606007642E-001 
1 .5237256040^0716-001 
1 .9 06968709 85O05E-001 
2.21641 112220419E-001 

2 • 43 1 1 476 7633960C-0OI 
2. 53390364 178345E-00* 
2 . 3 2376807632722E-00 1 
2. 4 198337O764297E-0OI 
2.24O6370531 1999E-001 
2.01994161 173123E-00I 
1. 769051 67599749E-O01 
1 . 37631 976499981 E-OOl 
1 .404853933876391E-001 
1 • 29003722671 106E-O01 
1.241 64 3762328*1 e-OOl 
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Claims 



5 1 . A method for the excitation by a radio frequency signal, in an NMR investigation apparatus, comprising 

electronic means for the treatment of the information in response to which the oscillations of the magnetic 
moments of nuclear spin contained in a predetermined volume to be examined submitted for the investigation 
are defined in amplitude and in phase for a set of given resonant frequencies, the said method being charac- 
terized by the following stages : 
10 — the said means are supplied with the Schroedinger equation relating to the probability amplitude of the 
spin wave function of these particles resonating at a given resonant frequency, the Hamiltonian of the said 
equation being constituted by a function dependent on the frequency and on the form of an excitation pulse 
applied to these spins, 

— the said Schroedinger equation is integrated with respect to time in the said means and the polynomial 
15 analytical expressions of the probability amplitudes are deduced at the end of the pulse as a function of 

the resonant frequency, said expressions being on the one hand of a degree equal to a number of elemen- 
tary periods whose sum of the durations is equal to the duration of the said pulse and on the other hand 
being modified by parameters constituted by coefficients representative of the effective excitation profile 
caused by the said pulse and submitted, during the same periods, by the spins at the different points in 
20 the said investigated volume, 

— using the said means to optimize the said coefficient in such a manner that the said polynomials make 
the best possible approach to a predetermined desired form, 

— deduction for each of the said elementary periods of the characteristics of an excitation signal leading 
to the said predetermined desired form, and 

25 — excitation of the spins of the said volume with this excitation signal. 

2. The method as claimed in claim 1, characterised in that in order to integrate an approximation using hard 
pulses is performed. 

3. The method as claimed in claim 1 or in claim 2, characterised in that for optimization the polynomial exp- 
ressions are investigated using minimum phase shaft, such expressions corresponding to a minimum power 

30 excitation. 

4. The method as claimed in claim 1 or claim 2, characterised in that excitation is performed using an exci- 
tation whose rate is similar to that of the representative figure 2 of the following table, the periods being given 
in a reduced time and the amplitudes being given in radians : 
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- 3. 904 9 1 06682062 1 C-002 
-6. *66*l 960134976C-O03 
-9. 2M 66596396 407E-O03 
1967746304O708E-OO3 
-6. 6J4 *26 23996 1 23E-003 
-6. 17591 456532324 E-003 
-7. 03961 632094723E<->003 
-3. 63477357 258307E-C03 
-3. 8e2339 1 101S924E-003 
- 1 . 7© 4 6 6*351 3643GE-003 

0. 234766360249696.004 
264 4 069O34366^E~OO3 

6 . 093630 7659 1 332E-OC3 
9 • 0321 6904390393C-403 
1 • 20553 1 4236 1 760C-002 
1 . 50 1 967 66 1 9 1 87 4 £-002 

1 . 7 65 1 74 1 33967 37C-0O2 

2. C43373463762S9£-002 
2. 27246 1 68283321E-002 
2. 436301 46956734E-002 
2.3071O287133642E-OO2 
1*. 6357 86* 1 203690E -O02 
2. 6343 1 355656 1 37 E -002 
2.5759 1 S 39192740C-002 

2. 41545403133147E-002 
2. 169569178037666-002 

1 . 637 1 3380632749E-002 
1.4J 9450621 09 303E-002 
9. 20457660665527E-r003 

3. 46765560646096E-OO3 
- 2. 92329 1 O0733393E-O03 
-9. e4?3390665»224E-^a 
-1. 715301 73784653E-O02 
-2. 4637153581 6664C-O02 
•3. 2 1 4422 1 94642 1 4C-002 
- 3. 936936 1 40 1 6602C-002 
-4 . e062 12 60029 1 67C-002 
-5. 19303174896276C-002 
-3. 66 36469 67486 1 OC-002 
-5. 9692228522574 9E-002 
-6. 1312131933939OE-O02 
•6.03961 140776231E-OO2 
-3. 744947391 12293E-002 
•5. 15620362699079E-002 
•4 « 26629999764 6 63E-002 
-3.05055661 136646E-O02 
-1 . 466 161 06373401E-002 

4. 48331 636317662E-003 

2. 773 1 4 1 67977S65E-002 
3. 5046 1 46234 17 1 SC-002 
e* 649 729O5264657C-O02 
t . 2193 1 56 1 4 22393E-O0 1 

1. 6OQ97O65O90379E-OO 1 

2. 0220306397 1342C-O0 1 
2.43839077421O71£*O0l 
2.8262393922193IC«-O0l 

3. 151 669792 18870C-001 
3. 36828780639445E-O01 
3. 44519 1 01 22O846E-O0 1 
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5. The method as claimed in claim 3, characterized in that for excitation an excitation is performed using 
an excitation whose rat is similar to that of the r presentativ figure 3 f the following able, the periods being 
given in a reduc d tim and the amplitud s being given in radians : 
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10 



Table 2 



-6. 4*6201 5321999 1 £-002 
-2. 22S06739796626E-002 
-2. 49869020943 1 08E-002 
-2. 68990157438880O002 
-2 • 7 7707 368058623E-002 
-2. 735947 i 191651 1 £-002 
-2. 54 4 794345579556*002 
-2. 1 863547248827 9E-002 
-1.64933463610OO9E-0O2 
-*.2975U66465l41E-0O3 
-3. 191 32207764 24 6C-004 

1 .030290774 06931 £-002 
2. 2324 950050628 1 E-002 

5.53908627237378E-002 
4. 90274 1906 1531 3E-002 
6.26431 148417868E-002 
7.5346637 1967290E-002 
8.69479228151243E-002 
9.601 552214151 16E-002 
1 .01921 802249669G-OO1 
1. 039 13 138439 167E-O01 
1 .01380187889778€-001 
9. 39 1 5552 1 524 1 7 8E-002 
2s ©• 
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6. 3739022 1 677407C-O02 
4. 1 37 90603273033E-002 

1.47486348149161E-002 
•1 . 54 8392344001 07E-O02 
-4 . G4 424 4 1 0 1 6884 1 £-002 
-Cr . 2985233 1684 852C-O02 
17625795521 197E-001 
<-1 .30471 46364 8678E-001 
-1. 79244 867442907£-00t 
-2. 01 439O13042493E-OO1 
-2. 147692056130606-001 
-Sr. 174347B4481 1 17E-001 
-2. 087 1 69 1 8900358E-001 
-1 ,8903O0?7066345£-001 
- 1 . 59 7389703706 1 6E-001 
- 1 . 22* 4 7 407 6864 39E-00 1 
-8 . 040576 1 297993 1 E-002 
-3. 436986O8033687E-O02 
1.34844833846243E-002 
6.16742832314041E-002 

1 . 00S826O6007&42C-OO1 

1.32372560401071E-O0t 
45 1 . 90696870985OO5E-O0 1 

2. 21 84 1 1 122204 1 9E-001 
2.431 14767835960C-001 
2. 33390364 178345C-001 
2.32576BO7632722C-0O1 
2.41 98 3 370764 297 £-00 1 
2.2408370531 1999E-001 
2.01994161 175123E-O01 
1 . 78903 1 6739 9749E-O0 1 
1 • 5765 1 976499981 E-00 1 
1.404 85393387081 E-OOl 
1 . 2900372267 1 1 06E-001 
1.241 643762328*1 E-001 
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